Optimal Scheduling for Fair Resource Allocation in 
Ad Hoc Networks with Elastic and Inelastic Traffic 
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Abstract — This paper studies the problem of congestion control 
and scheduling in ad hoc wireless networks that have to support 
a mixture of best-effort and real-time traffic. Optimization and 
stochastic network theory have been successful in designing archi- 
tectures for fair resource allocation to meet long-term throughput 
demands. However, to the best of our knowledge, strict packet 
delay deadlines were not considered in this framework previously. 
In this paper, we propose a model for incorporating the quality 
of service (QoS) requirements of packets with deadlines in the 
optimization framework. The solution to the problem results in 
a joint congestion control and scheduling algorithm which fairly 
allocates resources to meet the fairness objectives of both elastic 
and inelastic flows, and per-packet delay requirements of inelastic 
flows. 

I. Introduction 

As wireless networks become more prevalent, they will be 
expected to support a wide variety of services, including best- 
effort and real-time traffic. Such networks will have to serve 
flows that require quality of service requirements, such as 
minimum bandwidth and maximum delay constraints, while 
at the same time keeping the network queues stable for data 
traffic and guaranteeing throughput optimality. For the case of 
wireless networks with best-effort traffic only, optimization- 
based algorithms which naturally map into different layers of 
the protocol stack have been proposed in the last few years 
[l]-[6]; see [7] for a survey. However, these models do not 
take into account strict per-packet delay bounds. 

Scheduling packets with strict deadlines has been studied in 
[8]-[l 1], but all of these papers provide approximate solutions. 
The model that we study in this paper builds upon the recent 
work in [12]— [14] on admission control and scheduling for 
inelastic flows in collocated wireless networks, i.e., networks 
where all links interfere with each other. Among the many 
contributions in these papers is a key modeling innovation 
whereby the network is studied in frames, where a frame is 
a contiguous set of time-slots of fixed duration. Packets with 
deadlines are assumed to arrive at the beginning of a frame and 
have to be served by the end of the frame. In this paper, we 
explore this modeling paradigm further to study the design 
of resource allocation algorithms for ad hoc networks. The 
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frame-based model allows us to incorporate delay deadlines in 
the optimization framework for very general network models, 
and somewhat surprisingly, allows us to design a common 
framework for handling both elastic and inelastic flows. 
The main contributions of the paper are as follows: 

1) We present an optimization framework for resource 
allocation in a wireless network consisting of both best- 
effort flows and flows that generate traffic with per- 
packet delay constraints. The framework allows for very 
general interference, channel and arrival models. 

2) Using a dual decomposition approach, we derive an 
optimal scheduling and congestion control algorithm 
that fairly allocates resources and ensures that a required 
fraction of each inelastic flow's packets are delivered 
on time by appealing to connections between Lagrange 
multipliers, queues, and service deficits. The schedul- 
ing algorithm seamlessly integrates inelastic and elastic 
traffic into a unified max-weight scheduling framework, 
extending the well-known results in [15]. 

3) The convergence of the above algorithm in an appropri- 
ate stochastic sense is proved and it is also shown that 
the network is stable. 

II. Network Model 

The network is represented by a directed graph Q = (AT, £), 
where M is the set of nodes and L is the set of directional 
links such that for all n\, ri2 G N if (tlx, ri2) £ C then node 
ni can transmit to node ri2. The links are numbered 1 through 
\C\, and by abusing notation, we sometimes use I 6 C to mean 

Ze{l,2,...,|£|}. 

Traffic is assumed to be a mixture of elastic and inelastic 
flows, where an inelastic flow is one that has maximum per- 
packet delay requirements. In contrast, elastic flows do not 
have such requirements. 

Time is divided in slots, where a set of T consecutive time 
slots makes a frame. We assume that packet arrivals only occur 
at the beginning of a frame, and every inelastic packet has a 
deadline of T time slots. If a packet misses its deadline it is 
discarded, and it is required that the loss probability at link 
I G C due to deadline expiry must be no more than pi. For 
elastic traffic we associate a utility function Ui(xi) which is 
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a function of the mean elastic arrival rate per frame xi. We 
assume that J//(.) is a concave function. 

For a given frame, we denote by the vector a^ = (au)i^c 
the number of inelastic packet arrivals at every link, where an 
is a random variable with mean A; and variance a 2 a . We further 
assume that arrivals are independent between different frames 
and that Pr(a a = 0) > and Pr(au = 1) > 0. The last two 
assumptions are used to guarantee that the Markov chain we 
define later is both irreducible and aperiodic, although these 
can be replaced by other similar assumptions. Similarly, we 
define a e = (a e i)i e c to be the number of elastic packet arrivals 
at every link in a given frame. 

The channel state is assumed to be constant in a given 
frame, independent between different frames, and independent 
of arrivals. The vector c = (q);e£ denotes the number of 
packets link I can successfully transmit on a time slot in a 
given frame. 

Depending on the wireless technology used, we can have 
some channel feedback before or after a transmission occurs. 
If channel estimation is performed before transmitting, we 
can determine the optimal rate at which we can successfully 
transmit. Alternatively, feedback from the receiver after the 
transmission can be used to detect if a transmission is suc- 
cessful or not. In this paper we try to capture both scenarios 
in the following cases: 

1) Known channel state: It is assumed that q is a non- 
negative random variable with mean q and variance cr^, 
and we get to know the channel state at the beginning 
of the frame. 

2) Unknown channel state: It is assumed that c% is a 
Bernoulli random variable with mean q and we only 
get to know the channel state at the end of the frame. 

In the known channel state case where we do channel 
estimation to determine the optimal transmission rate, we 
can potentially send more than one packet in a time slot 
at higher rates. This is captured by the fact that we make 
no assumptions on the values q can take since it will be 
determined by the particular wireless technology used. In the 
case of unknown channel state we assume that we only get the 
binary feedback of acknowledgments, which is reflected in the 
Bernoulli assumption on q . In this case, and without any loss 
of generality, we assume only one packet can be transmitted 
per time slot per link. 

In the rest of the paper we will consider the known channel 
case first and then in Section [VI] we will highlight the 
differences in the analysis for the unknown channel case. 

III. Problem Formulation 

We first formulate the problem as a static optimization 
problem. Using decomposition theory, we will then obtain a 
dynamic solution to this problem and prove its stability using 
stochastic Lyapunov techniques. 

A feasible schedule s = (su,t> s e i,t) is such that 8u,u s e ;,t 
respectively denote the number of inelastic and elastic packets 
that can be scheduled for transmission at link I G C and time 
slot t 6 {1,2,..., T}; thus, su t t+s e i,t > means that link I is 



scheduled to transmit in time slot t of the frame. Furthermore, 
for any t, if s ih >t + s eh , t > and s ih >t + s e i 2 lt > then links 
1% and I2 can be scheduled to simultaneously transmit without 
interfering with each other. Assuming the inelastic arrivals and 
the channel state are given by a, and c respectively, we have 
the following constraints: 

T 

sn t t < an for alH G C and (1) 

t=i 

sa,t + s eUt < q for all / G C and t G {1, 2, . . . , T}. (2) 

We denote by S(ai,c) the set of all feasible schedules when 
the arrival state is a, and the channel state is c; thus, S(ai,c) 
captures any interference constraints we have on our network 
and satisfies (HJ and d2j. 

At the beginning of any frame we must choose a feasible 
schedule to serve all links and decide how many elastic packets 
are allowed to be injected in the network. Therefore, our goal 
is to find a function Pr(s\ai,c) which is the probability of 
using schedule s G S(ai,c) when the inelastic arrivals are 
given by a, and the channel state is c, subject to the constraint 
that the loss probability at link I G C due to deadline expiry 
cannot exceed pi. For elastic traffic, we want to select the 
vector a e such that we maximize the network utility while 
keeping the queues stable. 

To properly formulate the problem, let us first define 
Hi(ai,c) to be the expected number of inelastic packets served 
if the number of packet arrivals is given by a; and the channel 
state is c. Similarly, fi e (ai,c) denotes the expected number 
of elastic packets that can be served. Therefore, we have the 
following constraints: 

T 

fia(ai,c)< ^ 22 s u,tPr(s\ai,c) 
ses(ai,c) t=x 

T 

Hei{a u c)< ^ ^ s e i :t Pr(s\ai,c). 

seS(ai,c) t=l 

The expected service for mixed traffic at link I is then given 

by 

Mii = ^ ^ i± t i(ai,c)Pr{c)Pr{ai) 

di c 

li-ei d = ^ ^ Hei(ai,c)Pr(c)Pr(ai) 

a,i c 

and due to QoS requirements and capacity constraints, we 
require that 

fJ-u > A;(l - pi) and x t < /u e j. 

We will focus on maximizing the following objective for 
some given vector w £ R + ' : 

max y^ Uj{xi) + Wifiq (3) 

l*i(a,i,c),He(ai,c), 
[ii ./i e ,x,Pr(s\a,i ,c) 
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subject to 

T 

fiil(a>i,c) < J! S U,t Pr ( S \ a ii C ) for a11 1 G A a i' C 



s t=l 

T 



Hel(<H,c) < y^y^ Sei,tPr{s\a,i, c) for all Z e £, a,, c 

s i = l 

Mii = ^ ^ Mil ( a ii c)Pr(c)Pr(a,i) for all / € £ 
Mei = ^ y~l Vei(ai,c)Pr(c)Pr(ai) for all / € £ 

Mil > Aj(l - pi) for all Z e £ 
< xi < (j, e i for all Z G £ 
Pr{s\di, c) > for all s S 5(<Xi, c), aj, c 
y~^ Pr(s|ai, c) < 1 for all a^, c. 

s 

The vector w can be used to allocate additional bandwidth 
fairly to inelastic flows beyond what is required to meet 
their QoS needs. Other uses for w will be explored in the 
simulations section. We will assume that the arrivals and loss 
probability requirements are feasible and thus the optimization 
problem has a solution (x*,fi*). 

IV. Solution Using Dual Decomposition 

Using the definition of the dual function [16], we have that 

D(Si,S e ) = 

Ef Ui(xi) +WIHU - S et [xi - fiel] 
i x r ci \ i 

subject to 



lec 



Mil 



(ai,c) < X! s n,tP r ( s \ a ii c ) f° r ai l I £ A a i, c 



s t=l 
T 



fj, e i(ai,c) < X! X! s el,t-?M s l a i> c ) for all Z e £, a i; c 

s t = l 

fJ-il = / J 'ii( a i' c ) Pr i c ) Pr ( a i) for alH £ £ 

Mel = ^ fi e i(a i ,c)Pr(c)Pr(a i ) for all Z € £ 

> for all Z € £ 

Pr{s\a,i, c) > for all s S c), aj, c 
y~^ Pr(s|ai, c) < 1 for all aj, c. 

s 

Slater's condition [17] states that, since the objective is 
concave and the constraints are affine functions, the duality 



gap is zero and therefore D{8*,81) — Ui(x*) + wi/j,*^ 

lec 

where 

(S*,S*)e argmin D(6i,S e ). 

6 u >0,5 el >0 

We are interested in finding (x*,fj,*) but not the value 
D(S*,Sl), so if we rewrite the objective in the dual function 

as 

J2 Ui(xi) - SelXl 
lec 

+ S (Wl + Su)fJ,U + S e lH e l 

lec 

- Su^ii 1 - pi) 

lec 



max 

fj,i,^ e ,x ,Pr{s\di,c) 



we notice that the problem can be decomposed into the 
following subproblems: 



maxt/i(x;) - S e ixi 

x t >0 



and 



max S2(wi + 6u)[iu + 6 e in e i (4) 

/ Lij(oi,c),/i e (a ( ,c), £ 
fii,fi e ,Pr(s\ai,c) 

subject to 

T 

fMl(ai,c) < su,tPr(s\a,i, c) for all Z 6 £,di,c 

s i=l 
T 

Hei(ai,c) < s e i yt Pr(s\a,i, c) for all Z 6 £,ai,c 

s 1=1 

Mil = X! X! MiiK) c)Pr(c)Pr(ai) for all Z e £ 
Mei = 51 Mei(ai) c)Pr{c)Pr{a l ) for all Z e £ 

di c 

Pr(s\a,i,c) > for all s € S(a,i, c), cij, c 
y~^ Pr(s|ai, c) < 1 for all a*, c. 

s 

Furthermore, since we are interested in solving the problem 
for non-negative values of 5a and S e i, it must be the case 
that [i* and /i* are as large as the constraints allow, and since 
the upper bounds for n a (a,i,c) and ^(a^c) are expressed 
as a convex combination, and the objective function in 
is linear, the problem can be decomposed into the following 
subproblems for fixed ai and c: 



{T T 1 

t=l t=l J 



max 

ses 



This suggests the following iterative algorithm to find the 
solution to our optimization problem, where k is the step index 
and X max > max; e £ x\ is a fixed parameter: 

G argmax Ui(x t ) - 8 e i(k)xi 

0<x t <X max 
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s*(a;, c, k) G 



{T T ~j 

[w; + Sn (k)\ ^2 s u,t + $el (k) E s ei,t > 
t=l t=l J 

T 



a* c £= 1 



= E E E c ' fc)^(c)^r( ai ) 



flj c 1 



We update the Lagrange multipliers Si(k), 8 e (k) at every step 
according to the following equations: 



Sa(k + 1) = {^(fc) + e[Aj(l -pi) - 



and 



<5 ei (fc + 1) = {<*„,(*) + e[£?(fc) - ^(fc)]} + 



where e > is a fixed step-size parameter, and for any ael, 
a T = maxja, (J }. 

Making the change of variables ed(k) = Si(k) and eq(k) = 
8 e (k), we have that our iterative algorithm can be rewritten as 

x*(k) G argmax -Ui(x{) - qi(k)xi 
o<x,<x ma , e 



s*(a iy c, k) G 

f 1 - T T \ 

argmax^ I [~wi + di(k)} E s iM + E Se/ '* [ 



M^( fc ) = EEE s «.t( ai ' c ' fe ) Pr ( c )- Pr ( ai ) 

at c t—1 



W = E E E ~ S *lMi C ' 

C t — 1 



with update equations 

dt(k + 1) - [d,(fc) + A ; (l - ft) - M*;(fc)] + 

«(* + !) = &(*)+ 5? (*)-£rf(fc)] + - 

It should be noted that due to the change of variables d; (fc) 
can be interpreted as a queue that has A;(l — pi) arrivals and 
jj.*i{k) departures at step k; qi(k) can have a similar queue 
interpretation. The dual decomposition approach only provides 
an intuition behind the solution, but the real network has 
stochastic and dynamic arrivals and channel state conditions. 
In the next section, we present the complete solution which 
takes into account these dynamics and we also establish its 
convergence properties. 



V. Dynamic Algorithm and Its Convergence 
Analysis 

A. Scheduler and Congestion Controller 

To implement the algorithm online, we propose the follow- 
ing congestion control algorithm in frame k, where the queue 
length at link I is given by qi(k): 

x*(k) € argmax -Ui(x{) — qi(k)xi. (5) 

o<x,<x max e 

We need to convert this elastic arrival rate, which in general 
is a non-negative real number, into a non-negative integer 
indicating the number of elastic packets allowed to enter the 
network in a given frame. This conversion can be made in 
many different ways: we assume the elastic arrivals at link I, 
a e i(k), are a random variable with mean xf(k) and variance 
upper-bounded by a\, and are such that Pr{a e i{k) = 0) > 
and Pr(a el (k) = 1) > for all I G £ and all k. The last 
two assumptions are used to guarantee the Markov chain we 
define below is both irreducible and aperiodic, although these 
can be replaced by other similar assumptions. 

Letting the number of inelastic arrivals be denoted by ai(fc) 
and the channel state by c(fc), we propose the following 
scheduling algorithm: 



s*( ai (k),c{k),d(k),q(k)) G 



(6) 



arg max 

seS(oi(fe),c(ifc)) 



f 1 T T 

E \ ^~ wi + d i( k )\ E SiLt + E Sei ' t 

lec { 6 t=i t=i 



The vectors d(k) and q(k) are updated from frame to frame 
as follows: 

dt{k + 1) - [dt{k) + a a (k) - I^a^k), c(k),d(k),q(k))}+ 
q l (k+l) = [qi(k)+d el (k)-i: i (a i (k),c(k),d(k),q(k))} + , 
where 

T 

Z2(a,(fc), c(k),d(k),q(k)) = J2 ~S*l,Mk),c(k), d(k), q(k)) 



I* el ( ai (k),c(k),d(k),q(k)) 



1=1 

T 

E S *eiA a i( k ), c(k),d(k),q(k)) 
t=i 

and du(k) is a binomial random variable with parameters 
au(k) and 1 — pi- The quantity du(k) can be generated by 
the network as follows: upon each inelastic packet arrival, toss 
a coin with probability of heads equal to 1 — pi, and if the 
outcome is heads, add a one to the deficit counter. 

In our notation we make explicit the fact that for fixed e and 
w, the optimal scheduler (|6jl is a function of <2j(fc), c(k), d(k), 
and q(k). We interpret di(k) as a virtual queue that counts the 
deficit in service for link / to achieve a loss probability due to 
deadline expiry less than or equal to p\. This deficit queue was 
first used in the inelastic traffic context in [12] for the case of 
collocated networks; the connection to the dual decomposition 
approach now provides a Lagrange multiplier interpretation to 
it and allows the extension to general ad hoc networks. Note 
that qi{k) is just the queue size for elastic packets at link I. 
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B. Convergence Results 

For readability, we present the main results in this section, 
but the proofs are deferred to the appendixes. We start by 
noting that (d(k),q(k)) defines an irreducible and aperiodic 
Markov chain. To prove that our dynamic algorithm achieves 
the optimal solution to the static problem (O in some average 
sense and fulfills all links' requirements, we will first bound 
the expected drift of (d(k),q(k)) for a suitable Lyapunov 
function. 

Lemma 1: Consider the Lyapunov function V(d, q) = 
h^iectf + 9?- If /*« > Ml -pi) and M * ; > x* for all 
I e C, then 



E[V(d(k + l),q(k + l))\d(k) 
<Bi- B 2 Y,di- B 3 ^2qi 



d,q{k) 



V(d,q) 



1 



lec 



- E 



lec lec 
WiE[I^(ai(k),c(k),d, q)] 



-^mx^-UtixUk))] 



lec 



for some positive constants B\, B2, B 3 , any e > 0, where 
(x*,(jl*) is the solution to (0, x*(k) is the solution to (f5]), 
and I* (a,i(k) , c(k) , d, q) is obtained from the solution to (O. 

o 

It is important to note that since the last two terms in 
the right-hand side of the inequality can be upper-bounded, 
Lemma [T] implies that (d(k),q(k)) is positive recurrent since 
the expected drift is negative but for a finite set of values of 
(d(k),q(k)). As a direct consequence of this fact, we note 
that the total service deficit and queue length have a 0(l/e) 
bound. 

Corollary 1: If > A;(l - p{) and fi* el > x* for all / € 
C, then the total expected service deficit and network queue 
length is upper-bounded by 



lim sup E 

k — >oo 

for alH e C and 



.lec 



B 4 = 



< Ba 



1 



-Bk 



Bi 



and 



B s < 



mm{B 2l B 3 } 

J2iec max o<ii<x m „ 2\Ui(xi)\ + wiXi 
mm{B 2 ,B 3 } 



This also implies that the scheduling and congestion control 
algorithm fulfills all links' inelastic requirements. 

Corollary 2: If > Aj(l -pi) and [i* el > x ; * for all I e C, 
then the online algorithm fulfills all the inelastic constraints. 
That is: 



lim inf E 



1 K 

-^^( ai (fc), C (fc),d(fc), g (fc)) 



k=l 



>MI-Pl) 



for all I G C. o 
The above corollary simply states that the arrival rate into 
the deficit counter is less than or equal to the departure rate. 



This result is an obvious consequence of the stability of the 
deficit counters and so a formal proof is not provided here. 

Now we are ready to prove that our online algorithm is 
within O(e) of the optimal value. 

Theorem 1: For any e > 0, if > A; (1 — pi) and \i* el > x ; * 
for all I 6 C, then 



lim supi? 

K -too 



lec lec fc=i 

1 K 



K 

lec k=i 



< Be 



for some B > 0, where (x* , fx*) is the solution to (01, x*(k) is 
the solution to (O, and I* (oj(fc) , c(k) , d(k) , q(k)) is obtained 
from the solution to (0. o 
In conclusion, there is a trade-off in choosing the parameter 
e: smaller values will achieve a solution closer to the optimal, 
but at the same time the deficit in service at the links and the 
aggregate queue length increase. The statement and the proof 
of Theorem[T]follows the techniques in [3]. The result can also 
be derived, in a slightly different form, using the techniques in 
[4]. A closely related result can be obtained using the methods 
in [1]. 

VI. Unknown Channel State 

The analysis for the unknown channel case is similar to 
the one we presented for the known channel case, so in this 
section we will only highlight the differences. 

A feasible schedule s — (s«,{, s e i,t) is sucn that su,u 
Srfj respectively denote the number of inelastic and elastic 
packets that can be scheduled for transmission at link / G C 
and time t E {1,2, . . . ,T} without violating any interference 
constraints. Assuming the inelastic arrivals are given by Oj, 
and since we can only schedule at most one packet per link 
at every time slot, we have the following constraints: 



T 

E 

t=l 



su t < o-u for all Z € £ and 



(7) 



su,t + Sei.t < 1 for all I e C and t G {1, 2, . . . , T}. (8) 

We denote by S(a,i) the set of all feasible schedules for fixed 
arrivals, capturing any interference constraints we have on our 
network, and satisfying © and dS). 

Our goal now is to find a function Pr(s\a,i) which is the 
probability of using schedule s G S(ai) when the inelastic 
arrivals are given by <ij, subject to the constraint that the loss 
probability at link / G C due to deadline expiry cannot exceed 
pi. For elastic traffic, we still want to select the vector a e such 
that we maximize the total utility while keeping the queues 
stable. 

For a given distribution Pr{s\a{) we have that fiu(ai) is 
the expected number of attempted inelastic transmissions if 
arrivals are given by aj. Similarly, jU e i( a i) denotes the expected 
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number of times link I is scheduled to serve elastic packets in 
a given frame. As before, we have the following constraints: 

T 

M«( a i) - X! s H,t Pr ( s \ a i) 
s f.=l 

T 

^ei(ai) < y^y^ s e i. t Pr{s\a l ) 

s t = l 

When the (unknown) channel state is c, we have that 
cinuifli) is the expected number of successful inelastic trans- 
missions per frame at link I for fixed arrivals, while ci^, e i(ai) 
is the expected service to link I for inelastic arrivals. Thus, the 
expected service for mixed traffic at link I is given by 

Hil d = ^ ^ ciii i i(a l )Pr(c)Pr(a i ) 

c 

fJ-ei = y^y^c;^ e ;(ai)Pr(c)Pr(a i ). 

di c 

Simplifying both expressions we get 

m = y ^2,ciHii{ai)Pr{a i ) 

Mei = ^ cni e i{ai)Pr(ai). 

Due to service requirements and capacity constraints we 
need that 

fJ-u > A;(l -pi) and x t < \x ei . 

With the definitions and constraints stated above we can 
formulate the optimization problem in a similar way as in (0. 

The only difference with the known channel state case is the 
scheduling algorithm. Assuming inelastic arrivals are given by 
a,i(k) the scheduling algorithm is given by 

s*( ai {k),d(k),q(k)) G 

f 1 T T 1 

arg max ^ < [- wi + di (k)]c t ^ su,t + qi (k)ci ^ s e ;, t > ■ 
ses(a z (k)) leC [ e t=1 t=1 J 

The main difference in the scheduling algorithm compared 
to the known channel state case is that the network now uses 
the expected channel state in making scheduling decisions. 
Thus, the network needs to know or estimate cj as in [12]. 

Similar results can be proved for this algorithm using the 
techniques developed in Section IV-BI whereby one can show 
that the algorithm meets all the inelastic QoS constraints, the 
total expected service deficits and the queue lengths have a 
0(1/ e) bound, and the mean value of the objective is within 
O(e) of the optimal value. 




Fig. 1. Interference graph used in the simulations 



VII. Simulations 

The purpose of this simulation study is to understand 
how the parameter e and the link weights wi impact the 
performance of the algorithm, and how a greedy heuristic can 
be used to implement the optimal scheduler. We simulate a 
10-link network with an interference graph given by Fig. Q] 
where each node represents a link and each edge means that 
the two adjacent links cannot be scheduled simultaneously. For 
example, if link 1 is scheduled, then links 2, 4, and 7 cannot be 
activated. The required loss probability due to deadline expiry 
of inelastic packets is set to 0.1, the link arrivals are assumed 
to have a Bernoulli distribution with mean 0.6 packets/frame, 
and there are 3 time slots per frame. The channel for every link 
is assumed to have a Bernoulli distribution with mean 0.96, 
and we get to know the channel state at the beginning of the 
frame. We set Ui(xi) = log(x/) for all links. The simulation 
time was 10 6 frames. 

As can be noted from ©, the max-weight scheduler requires 
that we do an exhaustive search to find the optimal schedule at 
every frame. For large networks this can become a burden due 
to the large search space; thus we explore a greedy heuristic 
and check how close it is to the optimal solution: at any given 
time slot, the greedy scheduler orders all links according to 
their weights. The greedy scheduler adds one of the links with 
the largest weight to the schedule, then removes all links that 
interfere with this link from the graph, then schedules a link 
with the largest weight among the remaining links, and so on. 
This procedure continues until no more links can be scheduled. 

In Figs. |2] EJ and|4j we plot the expected values of the deficit 
counters and queues per link for various values of wi, and 
compare their evolution for both the scheduler with optimal 
decisions and the greedy scheduler. 

We see that as wi increases, the deficit counters become 
small. The upper bound in Corollary [T] only suggests that the 
sum of the deficit counters and queues is 0(1/ e). Thus, it 
is interesting to note that by changing wi, one can nearly 
eliminate the backlog in deficit for inelastic traffic while 
maintaining the same order of queue sizes. The reason for 
this can be understood by examining the scheduling algorithm 
©. Note that the algorithm gives priority to elastic traffic if 
queues are larger than counters. When wi is small compared 
to e, the effect of wi is negligible in the scheduling algorithm. 
On the other hand, when wi is 0(1), wi/e is 0(1/ e) which 
is comparable to the queue lengths and hence, the deficit does 
not have to be large to provide service to inelastic traffic under 
algorithm ©. 
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Fig. 4. Deficit size and queue length when tu; = 6 
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Fig. 3. Deficit size and queue length when ui; 



It must be noted that small deficit counters mean that there 
is a small backlog in providing acceptable service to inelastic 
arrivals. For the case of real-time traffic this is a desirable 
property, since we do not want to have large variations in the 
service provided that could affect the perceived quality. Thus, 
even if fair allocation of bandwidth beyond the minimum is not 
required for inelastic flows, choosing wi an order of magnitude 
larger than e is desirable to maintain small deficits. 

As can be noted, the greedy scheduler seems to give lower 
deficit values than the optimal scheduler for larger values of 
wi. We believe that the reason is that weights given to inelastic 
flows increase with increasing wi and therefore, the greedy 
scheduler picks them first. However, our optimality goal is 
given by (0 which is determined by the rates received by the 
various flows. The rates achieved by the two schedulers are 
quite close in the simulations, as seen in Figs.|3J|6l and|7l while 
keeping the dropping probabilities below the requirement, as 
shown in Figs. [S |U and [10] 



VIII. Conclusions 

In this paper we have presented an optimization framework 
for the problem of congestion control and scheduling of elastic 
and inelastic traffic in ad hoc wireless networks. The model 
was developed for general interference graphs, general arrivals 
and time-varying channels. Using a dual function approach 
we presented a decomposition of the problem into an online 
algorithm that is able to make optimal decisions while keeping 
the network stable and fulfilling the inelastic flow's QoS 
constraints. A key result is that, through the use of deficit 
counters, one can treat the scheduling problem for elastic and 
inelastic flows in a common framework. It is also interesting to 
note that the deficit counters introduced in [12]— [14] have the 
interpretation of Lagrange multipliers. Simulations corroborate 
our results and show the dependency of the performance of 
the algorithm on the auxiliary parameter e and its role into 
assigning resources to both elastic and inelastic traffic. We 
note that, in the unknown channel case, we do not consider 
channel feedback at the end of each time slot as in [12]— [14]. 
This will be addressed in future work. 
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Appendix A 
Proof of LemmaQ] 

To prove Lemma Q] we start by first proving two auxiliary 
lemmas and then stating a fact. 

Lemma 2: Given that at frame k we have the event d(k) = E 
d and q(k) = q, then 
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for some non-negative constant Bq, and where 



I*[(ai(k), c(k), d, q) is given by the solution to jSJ. 
Proof: 
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where (O follows from the definition of IZ(ai(k) , c(k) , d, q) 
and 
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Lemma 3: Given that at frame k we have the event d(k) 
d and q(k) = q, then 



E 



^{[ 9l +i el (i)-^(a,(fc),c(i),d, 9 )] + f 

<s 7 -i^(xn-^(^rw)] 

lec 

-^2<li{E [r el ( ai (k),c(k),d,q)] - x,*} 



E 
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Proof: 
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and ( fTTT > follows from the fact that is the optimal point 

of©. ■ 
Fact 1: The optimization in © can be performed over 
S(ai(k), c(k))cHi the convex hull of S(ai(k),c(k)); that is, 
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The reason for this comes from the fact that the objective 
function is linear and therefore there must be an optimal point 
s*(ai(k),c{k),d{k),q(k)))GS( ai {k),c(k)). ' ' o 
Proof of Lemma Q} For the purpose of this proof, we 
define the capacity region for fixed arrival and channel states 
a, and c as follows: 
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From Lemmas [2] and [3] we have: 
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Appendix B 
Proof of TheoremQ] 

From Lemma [TJ we know that 

e lec 

+ -^2wifi* l - wiE [J*;(ai(fc), c(k),d, g)] 
ze£ 

<s x - b 2 53 d ; - b 3 «l + V(d, q) 



lec lec 

- E [V(d(k + 1), q(k + l))|d(fc) = d, g(fe) = g] 
<Bi+F(d, g) 

- E [V(d(k + 1), g(fc + l))|d(fc) = d, g(fc) = g] 



since £>2 X)ze£ d i + ^3 Sze£ — 0- Taking expectations: 
e Lze£ 

+ E^/4 - wil* u (a t (k),c(k),d(k),q{k)) 
lec 

<B X - E [V(d(k + 1), g(fc + 1))] + £ [V(d(fc), «(*))] . 

Adding the terms for k = {1, . . . , K} and dividing by K 
we get: 



e 



WlV-u 



.lec 



E ^ E W(*0) + ^^K(fc), c(fc), d(fc), g(fc)) 
fc=i 

E[y(d(X + + , £[F(d(l),g(l))] 



K 

lec k=i 



<Bi- 
<B 1 + 



E[V(d{l),q{l))] 
K 



K 



(14) 



where ( TBI follows from the fact that the Lyapunov function 
V is non-negative. 

Assuming E [V (d(l) , q(l))] < 00 we get the following 
limit expression: 
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